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Abstract. We investigate interfacial structural and fluctuation effects occurring 
at continuous filling transitions in 3D wedge geometries. We show that 
fluctuation-induced wedge covariance relations that have been reported recently 
for 2D filling and wetting have mean-field or classical analogues that apply to 
higher-dimensional systems. Classical wedge covariance emerges from analysis of 
filling in shallow wedges based on a simple interfacial Hamiltonian model and is 
supported by detailed numerical investigations of filling within a more microscopic 
Landau-like density functional theory. Evidence is presented that classical wedge 
covariance is also obeyed for filling in more acute wedges in the asymptotic 
critical regime. For sufficiently short-ranged forces mean-field predictions for the 
filling critical exponents and covariance are destroyed by pseudo-one-dimensional 
interfacial fluctuations. In this filling fiuctuation regime we argue that the critical 
exponents describing the divergence of lengthscales are related to values of the 
interfacial wandering exponent deflned for planar interfaces in (bulk) two- 

dimensional (d = 2) and three-dimensional (d = 3) systems. For the interfacial 
height Iw ^ (9 — a)~^"' , with 8 the contact angle and a the wedge tilt angle, 
we find I3m = f(2)/2(l — C(3)). For pure systems (thermal disorder) we recover 
the known result (3^, = 1/4 predicted by interfacial Hamiltonian studies whilst 
for random-bond disorder we predict the universal critical exponent /3 ^ 0.59 
even in the presence of dispersion forces. We revisit the transfer matrix theory of 
three-dimensional filling based on an effective interfacial Hamiltonian model and 
discuss the interplay between breather, tilt and torsional interfacial fluctuations. 
We show that the coupling of the modes allows the problem to be mapped onto 
a quantum mechanical problem as conjectured by previous authors. The form of 
the interfacial height probability distribution function predicted by the transfer 
matrix approach is shown to be consistent with scaling and thermodynamic 
requirements for distances close to and far from the wedge bottom respectively. 
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1. Introduction 

Fluids adsorbed at micropatterned and sculpted solid surfaces may exhibit novel 
phase transitions and fluctuation effects compared to wetting behaviour at planar, 
heterogeneous walls A striking example of the influence of substrate geometry 

on fluid adsorption is provided by a simple wedge geometry j2| 13 El El E] ■ 
At two phase coexistence, a wedge-vapour interface is completely filled with liquid 
provided the contact angle 9 is less than the wedge tilt angle a. The phase transition 
from microscopic to macroscopic adsorption as — > a"*" is referred to as filling and 
may be first-order or continuous. The conditions for continuous wedge (and also conic) 
filling are less restrictive than for continuous (critical) wetting at planar walls 8 , 9 and 
give some hope that large-scale fluctuation effects associated with interfacial unbinding 
may be observable in the laboratory. Recent effective interfacial Hamiltonian studies 
have shown that fluctuation effects at continuous filling transitions exhibit a number of 
intriguing features. For three-dimensional wedge filling transitions critical singularities 
are believed to be far stronger than those characteristic of critical wetting transitions 
reflecting the anisotropy of soft-mode interfacial fluctuations induced by the wedge 
geometry. In particular for pure systems with sufficiently short-ranged forces the 
mean interfacial height l^ (measured above the wedge bottom) and roughness £,± are 
comparable and diverge with the same universal critical exponent = — 1/4 
IHl EI ■ These predictions are in very good agreement with Monte Carlo simulation 
studies of filling within more microscopic Ising and lattice polymer models |13[ [H] . 
For two-dimensional wedge filling on the other hand fluctuation effects are interesting 
for a different reason. Studies of fluctuation-dominated filling in both ordered (pure) 
[TH] and disordered (random-bond) JI] systems show that some observables, such as 
the mid-point height probability distribution function, show scaling properties which 
are identical with short-ranged critical wetting transitions. The only influence of 
the wedge geometry is to shift the effective value of the contact angle from Q to 
^ — a - a feature which has been referred to as wedge covariance jl2| . This "hidden 
symmetry" between wetting and filling appears to restrict the allowed values of the 
critical exponents at both 2D filling and wetting and leads to new some insights into 
the properties of critical wetting transitions. 

The present paper focuses on the structural and fluctuation properties of 3D 
wedge filling transitions. We begin with a discussion of wedge covariance and illustrate 
using a simple interfacial model of filling in shallow wedges that the fluctuation-induced 
covariance observed for 2D systems has a mean-field or classical precursor for wetting 
and filling in systems with short-ranged forces. Thus if Itt{0) denotes the mean- field 
result for the contact angle dependence of the critical wetting layer thickness, then in 
the wedge geometry the corresponding mean-field result for the mid-point height at 
bulk coexistence is 

U9,a)^lA9-a) (1) 

The predictions of the interfacial model are supported by a detailed numerical study of 
filling in a Landau-like density functional model that indicate this classical covariance 
is obeyed for both shallow and more acute wedges. The implications of classical wedge 
covariance for the structure of interfacial models of short-ranged filling and wetting 
are discussed. In the second part of our paper we turn our attention to fluctuation 
effects at 3D wedge filling and present a general scaling argument that relates the 
values of the critical exponents for 3D wedge filling to the wandering exponent of 
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planar, two-dimensional and three-dimensional fluid interfaces. Thus, for example the 
critical exponent (3^ is identified as 

(2) 
" 2(1 - C(3)) ^ ' 

where C,{d) is the standard wandering exponent for a planar-like interface in a d- 
dimensional bulk system. Our expressions recover the values quoted above for pure 
systems, corresponding to thermal fluctuations, and also allow us to discuss 3D flUing 
in disordered systems. We also discuss the cross-over from three-dimensional to two- 
dimensional wedge filling as one lowers the number of dimensions of translational 
invariance along the wedge. The change in the nature of fluctuation effects and 
critical singularities in this wedge compactification process highlights a remarkable 
numerical coincidence concerning the value of the wandering exponent for planar 
interfaces at the lower marginal dimension for wedge filling. In the final part of our 
study we focus on 3D filling in pure systems and address some problems that have 
been highlighted concerning the transfer-matrix analysis of a pseudo-one-dimensional 
effective Hamiltonian model. We argue that problems associated with the choice of 
the appropriate measure in the functional integral can be avoided if, in addition to 
breather-mode excitations of the interface, we also allow for tilt and torsional degrees 
of freedom. This more accurate formulation of the theory allows us to determine the 
universal scaling function associated with the density profile in the wedge. The form of 
this function is shown to have the correct short-distance and large-distance behaviour 
dictated by scaling theory and macroscopics respectively. 



2. Background theory 

2.1. Thermodynamics and critical exponents 

Consider a 3D wedge formed from the intersection of two smooth, planar walls that 
meet at angles a to the 2 = plane forming a wedge with opening angle tt — 2a. The 
parallel displacement vector in the z — Q plane is written x = (a;, y) with cartesians 
(a;, y) measuring distances across and along the wedge respectively. Thus the height 
of the wall above the z = plane is described by a height function = tana|x|. 

We suppose that the wedge is in contact with a three-dimensional bulk vapour at sub- 
critical temperature T and chemical potential tuned to bulk two-phase coexistence 
M — l^sat{T)^ . Provided the contact angle 0{T) of the sessile drop (defined for the 
planar wall-fluid interface) is less than ninety degrees the wedge preferentially adsorbs 
a volume of liquid near its bottom. Partial and complete filling refer to situations 
where the adsorption is microscopic and macroscopic respectively. The separatrix 
between partial and complete filling follows from simple thermodynamic arguments 
first discussed by Concus and Finn 4 some thirty years ago (see also Pomeau ^5^ and 
Hauge 0). Let V denote the accessible volume of the fiuid and L the area and 
length of the wedge respectively. The total grand potential Vl contains the macroscopic 
contributions 

n = -pV + (7u,vA + U,L (3) 

where a^v is the surface tension of the wall- vapour interface and is the excess wedge 
free energy. Now suppose we are at bulk coexistence and imagine that the wedge is 
filled to a height l^^, (see Figure^. Macroscopically the liquid- vapour interface must 
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be flat and hence must contain a thermodynamic contribution 

_ 2g;t,(cOSQ; - cos 9)lw , 
jw — ; (,4j 

sma 

where ai^, is the hquid-vapour tension and we have used Young's equation. It follows 
that for 9 < a the free energy can be lowered by completely filling the wedge. For 6 > a 
on the other hand the equilibrum value of 1^ must be finite and arises from a balance 
between thermodynamics and internal energy (intermolecular forces) or entropic 
(fluctuation) contributions. Thus by either increasing the tilt angle at constant T 
or increasing the temperature at constant angle (assuming the usual scenario where 
the contact angle decreases with temperature) one can induce a transition from partial 
to complete flUing at a temperature Tfm satisfying 

^{Tf^ll) = « (5) 

Filling transitions may be either first- or second-order corresponding to the 
discontinuous or continuous divergence of the mean interfacial height respectively. 
Importantly, the conditions for continuous filling are less restrictive than for continuous 
wetting transitions and wedges made from walls that exhibit first-order wetting 
(at some higher temperature T^et) may still exhibit continuous filling "SJ IHj. The 
most important diverging lengthscales which characterise the the transition are the 
mean mid-point height Iw, the perpendicular correlation length at the mid-point 
defined from the root-mean-square variance, and the correlation length describing 
fluctuations along the wedge. Correlations across the wedge are described by a 
lengthscale which is trivially related to the mean-height: ~ Iw cot a. At bulk 
coexistence critical exponents for the divergence of these length scales are identified 
according to 

L-^t-^- , U-t'"'- , iy-t'"' (6) 
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where the temperature- like scahng field t cx {Tfm ~T) cx (0 — a). The phenomenology 
is very similar for the case of 2D wedge filling except there is no analogue of the 
correlation length S^y. 

2.2. 3D wedge filling in pure systems 

A suitable starting point for the evaluation of the critical exponents is an effective 
interfacial Hamiltonian model based on a collective co-ordinate ^(x) measuring the 
height of the unbinding interface above the z = Q plane. At least for shallow tilt 
angles a tana, the appropriate effective Hamiltonian is a natural generalization of 
the standard capillary-wave model used to study planar wetting and we write [7] 

H[l]= JdKS^^{\/lf + W{l-a\x\)^ (7) 

where S is the stiffness coefficient of the unbinding interface which, for isotropic fluids, 
we can identify with the surface tension aiv . The model is not truly microscopic and is 
only valid for lengthscales much larger than the bulk correlation length. The binding 
potential W{1) describes the influence of intermolecular forces and decays at large 
distances as W{1) ^ —al~P. The Hamaker constant a is positive in the temperature 
region of interest whilst the value of exponent p depends on the specific range of the 
forces with p = 2, 3 for non-retarded and retarded van der Waals forces respectively. It 
is also possible to modify the model to include the effect of quenched disorder but for 
the moment we concentrate on pure systems with thermal fluctuations. Minimisation 
of the capillary-wave-like model Q determines the following mean-fleld values of the 
critical exponents for 3D filling |S| 

Pw^- , l'± = l , l'y = ^ + - (8) 

p 4 2 p 

Thus fiuctuation effects are extremely anisotropic at wedge filling with 3> S,x and 
are dominated by pseudo-one-dimensional local translations in the height of the filled 
region along the wedge jHHH- We refer to these as "breather-mode" excitations. From 
the values of the above critical exponents it is clear that mean-field theory is only valid 
for intermolecular potentials with p < 4 for which ^± <C Iw The filling fluctuation 
(FFL) regime corresponds to potentials with p > 4 and to study this Parry, Rascon 
and Wood introduced a pseudo-one-dimensional wedge Hamiltonian which accounts 
only for the breather- mode excitations jSHH: 

i?.N=/rf2/|^(^)' + nw} (9) 

where lo{y) = '(0, y) > is the local height of the interface above the wedge 
bottom and 1^ = {lo{y)). The model is considered valid only for small wave- vectors 
ky <^ fc™"^ ^ However because the fluctuations at fllling are strongly anisotropic 
the relevant scaling combination fyfc^"^ diverges in the scaling limit and the cut-off 
does not determine universal quantities. The most important feature of this effective 
model is the presence of a bending term resisting fluctuations along the wedge which 
is proportional to the local interfacial height. At bulk coexistence the wedge potential 
V{Iq) has the form 

Vilo)^-{e'-a')lo + Cll-P (10) 
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and ioT 9 > a has a minimum located at the mean-field value of l^^. Notice that 
the first term in the wedge binding potential is the small a. approximation to the 
thermodynamic term and is proportional to the linear, temperature-like scaling 
field t. The transfer-matrix analysis of the one-dimensional model will be considered 
later. Here we point out two related interpretations of the model which give the desired 
classification of the allowed critical singularities without need for explicit calculations. 
Both arguments will play a role in our later discussions. To begin we observe that the 
bending energy term is invariant under a renormalization group rescaling 

y^y' = l , io^i'o = ^ (11) 

with wedge wandering exponent C«j = 1/3- Under this transformation the linear 
scaling field t and wedge Hamaker constant C rescale to t' and C with 

t' = bh , C' = b^C (12) 

Thus as anticipated the intermolecular forces are relevant and irrelevant for p < A and 
p > 4 respectively. The temperature-like scaling field is always relevant and from its 
rescaling and the value of Ciu we can identify the values of the exponents in the FFL 
regime as 

113 

^w = -^ , ^-1 = 4 ' ^!/ = 4 (1^) 

Equivalently we may say that in the FFL the diverging lengths scale according to 

~ a ~ ^1 (14) 

which shows the dramatic influence the wedge geometry has on the fluid interfacial 
properties near the filling phase boundary. Closely related to this simple rescaling 
analysis is the construction of an effective potential. Heuristically we may imagine 
that the equilibrium interfacial height 1^, follows from minimization of an effective 
wedge potential which takes into account the contribution arising from the bending 
free energy. Noting the scaling relation between 1^, ^± and £,y together with the first 
term in the wedge Hamiltonian it is natural to suppose that the effective potential has 
the form 

Veff{lo)=Vilo)+Dl-^- (15) 

where =2/^^ — 3 and the constant D oc T,. Thus for the present case of thermal 
disorder = 3 and the critical behaviour falls into two regimes for p < 4 and p > 4 
depending on whether intermolecular forces or the bending energy is the leading order 
correction to the thermodynamic term. 

It is straightforward to generalise these ideas to a d-dimensional wedge which 
has translational invariance in = d — 2 dimensions along the wedge. As we shall 
see consideration of the "compactifaction" from a three- to two-dimensional wedge 
highlights some interesting properties of the wandering exponent for planar interfaces. 
Breather-mode excitations do not lead to large-scale interfacial roughness for d > 4 
and mean-field theory is valid. For d < 4 we find two possible fluctuation regimes 
corresponding to mean-field-like and fiuctuation-dominated behaviour. The FFL 
occurs for sufficiently short-ranged potentials with p>2(l/<^^ — 1) with a generalised 
wedge wandering exponent 

= ^ (16) 
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valid for 2 < d < 4. The critical exponent for the divergence of the interfacial height 
in the FFL regime follows as 

and recovers the 1/4 power law for the three-dimensional wedge. Intriguingly the result 
p7|l has the correct two-dimensional limit = 1 for filling with thermal fluctuations. 
We shall return to this later but first focus on the 2D limit in more detail. 

2.3. Fluctuation-induced wedge covariance in 2D 

Interfacial Hamiltonian studies ^1 ^2 reveal that the fluctuation regimes for 
2D filling are different to those for 2D critical wetting. Nevertheless there appears 
to be a profound connection between the fluctuation-dominated regimes for each 
transition. Recall that critical wetting refers to the continuous unbinding of the 
liquid-vapour interface (say) from a planar wall as the temperature is increased to 
the wetting temperature T^et at bulk coexistence. The vanishing of the contact angle 
9 as T —^ T^^et is accompanied by the divergence of the mean interfacial thickness 
?7r, roughness and transverse correlation length ^||. More detailed information 
concerning the local height fluctuations is contained in the one-point interfacial 
height probability distribution function (PDF) P^(/;0) where, in preparation for our 
discussion of wedge covariance, we have written the PDF in terms of the contact angle. 
Wetting transitions are characterised by highly anisotropic fluctuations arising from 
the random- walk or capillary-wave- like motion of the interface which are quantified by 
the scaling relation ~ C^^'''' with (^(d) the wandering exponent for a free interface 
(see later). At bulk coexistence the critical exponents for critical wetting describing 
the asymptotic divergence of these lengthscales are identified from 

I. ~ (T^et - r)-'^= , e± - {T^et - T)-"^ , ClI - {Tn,et ~ (18) 

Generically, fluctuation effects at critical wetting fall into three regimes depending 
by the range of the intermolecular forces |15[ I16L 1171 1181 [TU| . These are mean- field 
(MF), weak-fluctuation (WFL) and strong- fluctuation (SFL) scaling regimes. The 
SFL regime represents the universality class of critical wetting with sufficiently short- 
ranged forces p > 2(1/^(2) — 1) and is characterised by large fluctuations /tt ~ C-L and 
universal critical exponents. In two dimensions rather general random-walk-based 
arguments predict that the critical exponents in the SFL regime are explicitly related 
to the value of the wandering exponent. In particular for the interfacial height it is 
believed that [TTIIT^ 

3 = (19) 
1-C(2) ^ ' 

which is certainly obeyed for pure (C(2) = 1/2) and random-bond systems (C(2) = 
2/3). 

For 2D filling on the other hand interfacial Hamiltonian studies show there are 
only two fluctuation classes corresponding to mean-field and fluctuation-dominated 
regimes for which ^ £^± and ^ respectively |ll)llll[[T^ . In general the critical 
exponents describing mean-field filling and mean-field critical wetting transitions are 
unrelated. Also the filling fluctuation (FFL) regime, representing the universality class 
of systems with short-ranged forces, is broader than the SFL regime for critical wetting 
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and occurs for p > (1/C(2) — 1). Simple scaling arguments indicate that within the 
FFL regime the interfacial height diverges with a critical exponent 

0. - J% (20) 

which is in agreement with model calculations in pure |1U| and impure systems. 
Wedge covariance |12| refers to the "empirical" observation emerging from effective 
interfacial Hamiltonian and Ising model studies that the equality of the critical 
exponents also extends to the full scaling properties of the respective one-point PDF's. 
In brief, wedge covariance states that for systems with short-ranged forces, at bulk 
coexistence and in the scaling limit, the interfacial height PDF at the wedge mid-point 
Pw{l;6,a) is identical to the analogous PDF for strong- fluctuation regime critical 
wetting at a planar wall-fluid interface but with an effective shifted contact angle 
6 ^ 9 — a with a the wedge tilt angle. Thus if Pt^{1\9) denotes the critical wetting 
interfacial height PDF written in terms of the contact angle, wedge covariance implies 

P,S;0,a)=P^{l;e-a) (21) 

which is valid in the critical regime 9 — a —t 0. From this it follows that the equilibrium 
mid-point height near a fluctuation-dominated filling transition satisfies 

l^{9,a) = l,{9-a) (22) 

where, in an obvious notation, iTrifi) denotes the SFL regime critical wetting fllm 
thickness expressed in terms of the contact angle. Interestingly, the covariance (|21|l 
necessarily leads to the exponent identifications (|19I2U|) whilst a third relation between 
the FFL wedge free energy and SFL point tension t{9) leads to the new result for 
wetting [HI 

^.W = -^ (23) 

This expression leads directly to the identification of the critical singularity associated 
with the point tension consistent with exact Ising model calculations |2U| and a more 
general conjecture due to Indekeu and Robledo We emphasise that the above 

covariance relations, originally noted from interfacial Hamiltonian studies of filling 
in shallow wedges (in pure and impure systems), appear to be rather robust. For 
pure systems they are also obeyed by a drumhead interfacial model of filling in more 
acute wedges |22| and, most importantly, are consistent with exact |23| and numerical 
studies H3] of filling in the square lattice Ising model with tilt angle a = tt/A. 



3. Classical Wedge Covariance for filling in shallow wedges 

The standard fiuctuation theory of wetting at planar wall-fluid interfaces is based on 
analysis of effective interfacial Hamiltonian models which describe the fluctuations 
of a collective co-ordinate Z(x) measuring the local height of the interface from the 
wall. This coarse-grained description is valid at lengthscales much bigger than the 
bulk correlation length and, in its simplest form, is written 

H.[l]= j d^^{Vlf + W{l)^ (24) 

where, as described earlier, E is the stiffness coefficient of the unbinding interface and 
W{1) is the binding potential which, in general, accounts for the direct influence of 
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intermolecular forces. We shall focus on isotropic fluid interfaces for which S can be 
identified with the surface tension aiy . 

For systems with short-ranged forces the binding potential describes not so much 
the range of the intermolecular potentials but the decay of perturbations in the local 
microscopic order parameter (density, magnetization) and is usually taken to have the 
form, at bulk coexistence, 

W{1) = -ae'""^ + be-^''^ (25) 

where k is the inverse bulk correlation length of the wetting phase. The temperature 
dependence of the coefficients a, b is crucial at mean-field level. The leading order term 
vanishes at the mean-field wetting temperature so that a cx (T^ef — T) whilst b must 
be positive finite at T^et to ensure stability and is usually taken to be constant. The 
structure of the binding potential can be inferred by comparison with more microscopic 
Landau-like mean-field theories. More formally it can derived from a constrained 
fluctuation sum as in the approach of Fisher and Jin which also leads to the presence 
of a position-dependent stiffness-coefficient |25l I26| . 

At mean-field level the equilibrium thickness It^ of the interface follows from 
minimization of W{1) whilst the contact angle 9 can be identified from S6'^/2 = 
—W{lTr). Thus for systems with short-ranged forces the mean-field results are 

Kl^=\n2b/a , a ^ V2U)9 (26) 

corresponding to critical exponents f3s — 0(ln) and = 0. Written in terms of the 
contact angle the equilibrium film thickness of the critical wetting layer is therefore 

>il^ie) = -\n^e (27) 

which will be our intial point of reference for the wedge calculation. 

For shallow wedges corresponding to tilt angles a ~ tana the natural 
generalization of the planar model is the effective Hamiltonian jT)) which simply 
assumes the interfacial interaction with the wall is controlled by the local relative 
height variable / = ^ — a|x|. Note that the mean- field analysis is essentially 
independent of dimension since it assumes translational invariance along the wedge. 
The Hamiltonian is minimised subject to the appropriate boundary condition I ?^ 
as \x\ — > oo and yields the Euler-Lagrange equation 

i:i"{x) = W'{l-a\x\) (28) 

It is convenient to write this in terms of the local relative height I — I — a\x\ which, 
making analogy with classical mechanics, exploits the local Galilean invariance of the 
Euler-Lagrange equation. On integration one obtains the "energy" equation 

^i\nx)\~af^AWil-a\x\) (29) 

where AW{1) — W{1) — W{lTr). Thus at the wedge mid-point one finds the simple 
expression for the local height of the filling film 

= AW{L) (30) 

As pointed out by Rejmer et al this equation has a very elegant graphical 
interpretation which demonstrates that for quite arbitrary choices of binding potential 
the wedge undergoes a filling transition when 9{T) = a, in precise accord with 
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thermodynamic arguments. Now consider the specific case of short-ranged forces with 
the binding potential H25|) . Substitution determines the mid-point height as 

Kl^i9,a)^-\n^{9-a) (31) 

This not only identifies the logarithmic divergence of the film thickness at the filling 
transition, /?/ = 0(ln), but also reveals that the mean- field theory shows a "classical" 
analogue of wedge covariance observed in the 2D calculations: 

L{9,a)=l^i9-a) (32) 

In other words, for systems with short-ranged forces, the influence of the wedge 
geometry as manifest in the mid-point height is to shift the effective value of the 
contact angle. A similar property extends to the whole equilibrium profile which 
satisfies 

l{x) - a\x\ = l^{9 - a + \l'{x)\) (33) 

and smoothly interpolates from lTr{9 — a) to IjriO) as |a;| increases. Profile covariance 
indicates that the height dependence of the meniscus contains information about the 
contact angle dependence of the planar wetting film ?^(0) and vice versa. 

Classical wedge covariance is also manifest in the Gaussian fluctuations about the 
mean-field solutions. This is most easily interpretated for the 2D wedge. Consider for 
example the mean-field expression for the connected height-height correlation function 

suo,x) = im - mmm ~ m))) m 

and note that S'uj(0, 0) = identifies the mid-point roughness. The correlation 
function can be easily obtained from solution to the appropriate Ornstein-Zernike-like 
equation yielding 

where, for convenience, we have set ksT — 1. For systems with short-ranged forces 
the explicit result for the mean-field roughness at the wedge mid-point satisfies 

2kl,[0 — a) 

This can be compared with the result pertinent to critical wetting at the planar wall- 
fluid interface which can be obtained by simply setting a = 0, 

implying the classical wedge covariance relation 

C±{0,a)^Ui0-a) (38) 
for Gaussian fluctuations about the mean-field solution. 

For higher-dimensional wedges the actual mid-point roughness ^± is no longer 
covariant because of the breather-mode excitations along the wedge. For such 
dimensions the quantity analogous to the roughness which shows covariance is simply 
the zeroth moment of the the height-height correlation function 

S^iO,x) = / dy^^am - (?(0)))([(x) - ([(x)))) (39) 
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where t/12 is the relative separation of the two points along the wedge. This satisfies 
the same expression quoted above |S] and is therefore covariant. 

We emphasise that classical covariance is not a general feature of the mean-field 
theory of filling and critical wetting. Within a more general description based on a 
binding potential 

Wil) = ~al-P + br" (40) 

the filling and critical wetting exponents are distinct: 

Pw = - , A = ^ (41) 

p q-p 

This automatically rules out the possibility of (classical) wedge covariance for binding 
potentials of the form H40|l . Nevertheless the results presented above showing wedge 
covariance for short-ranged forces do generalise to the class of potentials 

W{1) = -auj{l) + buj{lY (42) 

where ix'(^) corresponds to an arbitrary choice of monotonically decaying function. 
Choosing w = e""' we obtain the usual short-ranged binding potential whilst setting 
u> = l~P one obtains a binding potential describing a particular type of multicritical 
wetting transition with q ~ 2p. Repeating the analysis above we obtain for the planar 
critical wetting and mid-point wedge filling interfacial heights 

'^a^)-/!^ ' (43) 

which immediately implies the wedge-covariant relation (|32|l . A little more algebra 
shows that the roughnesses also obey covariance implying that the mean- field PDF's 
are identical provided we map 9 9 — a. We refer to the class of potentials (|42|l 
as classical wedge covariant binding potentials. Of course, only for the short-ranged 
case u! = e~'^' do we anticipate that such potentials have any physical significance. 
Nevertheless they do point to an important feature of the present analysis. It is easy 
to demonstrate that the wedge covariant potentials 1)42(1 all describe planar critical 
wetting transitions with vanishing specific heat exponent = 0. This mirrors 
precisely the situation for the fluctuation-induced non-classical covariance since for 
C > 1/2 the value of the specific heat exponent in the SFL regime is as = 0. This 
covers both cases, C = 1/2 and ( = 2/3 in 2D, where non-classical covariance is 
known to occur and it is tempting to speculate that the vanishing of the specific 
heat exponent either at mean-field level or beyond plays a key role for classical and 
non-classical covariance respectively. 

The prediction of classical wedge covariance shows that this hidden relation 
between filling and wetting for special types of forces is not necessarily a fiuctuation- 
induced phenomenon. Having said that we feel care should be taken in saying that the 
classical covariance is the underlying origin of non-classical covariance. Whilst they 
are certainly related it is far from obvious why the covariant relation H32|l obeyed for 
binding potentials describing mean-field critical wetting with as = should remain 
unaltered in the presence of large-scale fluctuation effects which renormalise the value 
of the exponent f3s whilst leaving as unchanged. In view of this we treat classical 
covariance as a prediction which should be tested in more microscopic models of 
wetting and flUing. Since the binding potential l|25|) is believed to describe wetting in 
systems with short-ranged forces we turn attention to numerical studies of filling based 
on a Landau-like density functional model. This will also allow to test any limitations 
arising from the shallow wedge approximation implicit in the interfacial model. 
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4. Wedge Filling within Landau Theory 

For our Landau theory study we resort to a magnetic terminology rather than 
the fluids-based one considered earlier. At mean-field level the equilibrium order 
parameter m(r) is translationally invariant along the wedge so we can restrict ourselves 
to magnetization profiles in a two-dimensional space r — (x, z) with x the co-ordinate 
across the wedge. The free-energy frmctional for the infinite wedge that we wish to 
minimise is 

F\m\— f dr I -(S/m)'^ ~ -m"^ + —m'^ — hm 

where the volume of integration is restricted to z > tana|a;| for every x. The 
parameter t measures the deviation from the bulk critical temperature (which is 
always finite) whilst u > for stability. The bulk field h = 0~ so that the bulk 
magnetization is negative. The temperature dependence of the equilibrium profiles 
can be eliminated by measuring the magnetization in units of the bulk spontaneous 
magnetization mo(i) = \ftju and need not be specified further. Rather than use a 
local surface field and enhancement parameter we use fixed boundary conditions which 
set the surface magnetization to a positive value m(a;, tan a|a;|) = m\ for all x. This 
is equivalent to the infinite surface enhancement limit in the model of Nakanishi and 
Fisher EH and ensures that the wetting transition pertinent to the planar wall-down 
spin interface is always second-order. In the planar limit a = the model can be solved 
analytically and exhibits a critical wetting transition when the surface magnetization 

m^"* = mo (45) 

which allows us to induce wetting (and filling in the wedge) by either increasing m\ 
at fixed t or vary t at fixed m\. We have chosen to vary m\ sX fixed temperature since 
this keeps the bulk magnetization and correlation length ^j, = 1/k — (2^)"^^/^ fixed. 
The contact angle within this model can be calculated analytically 

cos6'=— 1-— 4j 46 

so that near the wetting transition cx (toq — m\)lmQ where (mo — mi)/mQ 
may be regarded as the temperature-like linear scaling field. Thus in the wedge 
geometry it is straightfoward to convert the numerically-determined value of the 
surface magnetization at filling phase boundary m{*" into a contact angle 6. We 
also remark that near the filling transition the scaling field 6 — a is equivalent to 
(m{*" — mi)/mo. 

We have numerically minimised a discretized version of the continuum Landau 
free-energy functional in a number of finite-size geometries with boundary conditions 
chosen to mimic the bulk as closely as possible. We considered a double wedge 
geometry with opposing surface magnetizations on the top and bottom wedges and 
also a capped wedge geometry (see Figure |2l with both free and fixed boundary 
conditions at z = Li. The results for the phase boundary and logarithmic growth of 
the filling layer are the same in all of these geometries. After testing our algorithm 
on a number of standard one- and two-dimensional problems including the Laplace 
and Poisson equations we first applied the minimization scheme to the Landau free- 
energy functional with a = 7r/4. We found that, as the system was made larger and 
the discretization finer, our results stabilised and the phase boundary could be made 
arbitrarily close to the theoretical prediction. The smallest system with the roughest 
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Figure 2. Typical capped wedge geometry used for the Landau numerical 
calculations. The magnetisation has a fixed value mi at the wedge boundaries 
and the bulk value mo at z = L\. Here, a = 45° and Li = L2 ~ SO^;,. Two 
solutions corresponding to either side of a filling transition are shown: the lower 
interface for rrti/mo = 0.5 and the upper interface for mi/rriQ = 0.55. 

discretization that produced a clear logarithmic growth and a transition value of 9 
within 0.1 degrees of a was Li = L2 ~ 30^b and a distance between the grid points of 
approximately Taking this as a starting point we increased the number of points 
on the grid and reduced the grid spacing in the vertical direction, producing a more 
open wedge of the same depth, with the same discretization in the horizontal direction 
and a finer discretization in the vertical. The number of points and the vertical scale 
were adjusted to produce the various wedge angles. A similar procedure was used 
to produce acute wedges, but here the horizontal distance L2 was kept constant and 
the wedge made gradually deeper. This approach allowed us study filling in wedges 
with a fairly wide range of tilt angles, up to a maximum of about 70 degrees, in a 
reasonable amount of CPU time. The effect of the variable discretization can be seen 
in Figure|31 For the most open wedges, 9 and a agree to four significant figures as the 
fine discretization in the direction of the film growth, whilst for the very acute wedges 
6 and a may differ by as much as 0.5 degrees. 

The results of our study are presented in Figure 13 and Figure ^ In Figure |2I we 
present data for the numerically determined phase boundary for tilt angles between 15 
and 75 degrees which is in excellent agreement with the theoretical prediction 9 = a. 
We have also studied the meniscus profile defined as the surface of isomagnetization 
m = (see Figure 01). For wedges with tilt angles a < 7r/4 there is clear evidence for 
the logarithmic growth of the mid-point filling height 

k/^ = Aln(mf"-mi) + C (47) 

with a universal, angle independent constant A = 1.02 ± 0.04 consistent with the 
shallow wedge limit A = 1. The constant C is non-universal as is the size of the 
asymptotic critical regime which decreases with increasing wedge angle. Nevertheless 
even for the most acute angles we consider there is no evidence that the amplitude A 
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Figure 3. Plot of the contact angle at the filling transition, 9, against the 
wedge angle, a. The error bars on 9 lie within the circles. The continuous line 
corresponds to the theoretical prediction Equation O . 




Figure 4. Plot of the reduced midpoint interface height, Z^/Cft, against the 
reduced surface magnetisation, mi /mo, for a range of a between approximately 
15° (bottom) and 75° (top). The dashed line corresponds to a = 45°. 
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Figure 5. Substrate-interface interaction paths (daslied lines) for the wedge 
filling transition in the model of Rejmer et al 



differs from unity. More precisely for all wedge angles where there is clear logarithmic 
growth of the filling layer the amplitude of the logarithmic growth is consistent with 
A = 1. In other words, even away from the shallow wedge limit the film diverges as 

K/„(6l,a) - -ln(6'- a) (48) 

and behaves, in the critical regime, precisely like /^(6'— a). This observation is certainly 
aesthetically appealing. In 2D, all results point to non-classical covariance for both 
shallow and more acute wedges as indicated by exact Ising studies for the right-angle 
corner and also a drumhead model calculation with arbitrary wedge angle a. The 
Landau theory numerics, coupled with the interfacial model studies for the shallow 
wedge, show that this is also the case for classical covariance. 

Given that the interfacial model Q is only valid for shallow wedges its is natural 
to inquire whether one can explain the Landau theory results using a generalised 
interfacial model valid for more acute wedges. We finish this section by showing 
that a very reasonable, and for many purposes satisfactory, drumhead-like interfacial 
model proposed by Rejmer, Dietrich and Napiorkowski (RDN) does not fulfil this 
requirement. The model is a straightforward generalization of the shallow wedge 
Hamiltonian but is amended in two ways. First the square gradient term multiplying 
the tension is replaced by the correct expression for the total area of the liquid- vapour 
interface. This is familiar from older drumhead-like capillary-wave models and ensures 
that the model recovers the correct filling phase boundary 9 = a. Secondly they 
propose replacing the relative vertical height from the wall to interface, appearing in 
the generalised binding potential contribution, by the normal distance to the closest 
wall (see Figure I^J. Thus the RDN model for the 2D/3D wedge is given by [7] 

+ sec aW {sec a [I - i/'(x])| (49) 

where, as earlier, I — /(x) denotes the vertical height relative to the x = plane 
and Tpi^) = tan a I a; I . This assumption is certainly valid far from the wedge bottom 
where the interface and wall are parallel. One potentially unsatisfactory feature of it 
is that points on the wall for which for distances \x\ <lw sin a cos a do not contribute 
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to the binding potential i.e. the interface is bhnd to the wedge bottom. The mean- 
field equation for the mid-point height, obtained by minimising the Hamiltonian and 
integrating once, is 

cr;t,(cos6' — cosa) = iy(^^ cosa) (50) 

which is the analogue of the shallow wedge expression H28|l . In this way it 
is straightforward to see that the wedge filling transition occurs at the correct 
thermodynamic phase boundary 9 — a and that the critical exponents for continuous 
filling are independent of the wedge angle. Now consider the case of short-ranged 
forces modelled by the binding potential (|25|l . According to the above drumhead 
model the asymptotic divergence of the mid-point height is given by 

kIw ^ — sec ahi{9 — a) (51) 

and contains a geometrical amplitude factor sec a which is not consistent with the 
Landau theory numerics. In other words the RDN model does not give the correct, 
wedge-covariant result H48|) . 

The failure of the RDN model to account for classical wedge covariance in more 
acute wedges is a surprising result. The model is certainly plausible and would appear 
to be the simplest possible generalization of the shallow wedge model. Indeed when 
considered beyond mean- field theory using transfer matrix methods, the RDN model 
does predict non-classical wedge covariance for 2D systems with strictly short-ranged 
forces! We emphasise that there is no mystery here; non-classical wedge covariance is 
a fluctuation-induced phenomenon for which the precise form of the binding potential 
(determining the way in which short-range forces are modelled) is irrelevant. In the 
presence of large fluctuations with 1^ ^ £,± all points on the interface can feel the 
influence of the wall through collisions with it. Thus in the 2D FFL regime it is 
sufficient to use a binding potential with a simple square well shape representing pure 
contact forces. This contrasts with classical wedge covariance which reflects the precise 
form of the underlying interfacial Hamiltonian model. At mean-field level there are no 
fluctuation effects that take the interface to the wall and even for short-ranged forces 
the manner in which one models the large distance exponential tail of the binding 
potential is crucial. Clearly the assumption that the wall-interface interaction occurs 
via the normal distance to points on the closest wall is incorrect. The absence of a 
sec a prefactor in the Landau numerics H48() indicate that the correct measure of this is 
more akin to an effective local, vertical interaction similar to the shallow wedge model. 
This observation may well have ramifications for the construction of interfacial models 
of 3D wetting and filling in systems with short-ranged forces. 

5. 3D wedge filling and the interfacial wandering exponent 
5.1. Interfacial wandering and criticality at complete wetting 

In the treatment of fluctuation effects at 3D wedge filling presented in the introduction 
there appears to be little connection between the values of the critical exponents 
and the interfacial wandering exponent ^(d) defined for a planar-like fluid interface. 
The value of the wandering exponent (^{d) is crucial for discussions of dcpinning as 
well as complete and critical wetting transitions. Interfacial fluctuation effects at 
such transitions are isotropic in the d — 1 dimensions parallel to the interface and 
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the wandering exponent describes the scahng relation between the roughness and 
transverse or paraUel correlation length ^|| |17| : 

U - 4^'^ (52) 

In general the value of C(d) depends on the dimensionality of space and the qualitative 
type of disorder. For pure systems, thermal fluctuations lead to interfacial roughness 
for d < 3 and we have the well-known result jl7| 

Cid) = ^ (53) 

where for the marginal case, the identification ^(3) = corresponds to ^ ^In^n . 
The lower critical dimension for phase separation is = 1 at which C.{d) = 1. 
The wandering exponent is altered in the presence of quenched impurities and it is 
usual to distinguish between random- field and random-bond disorder. Random flelds 
induce interfacial roughness for d < 5 and general scaling arguments lead to the result 
^(d) = (5 — (i)/3 |17[I28[EU| implying a lower critical dimension at — 2. For random 
bonds on the other hand the explicit dimension dependence of C(c?) is not known. Near 
the upper marginal dimension, d = 5, approximate functional renormalization group 
calculations of an interfacial model lead to the linear relation ^| ISU] 

C(d) « 0.2083(5 - d) (54) 

This is believed to be accurate even far from the marginal dimension since for d = 2 it is 
close to the known exact result ^(2) = 2/3 TP. In three dimensions numerical results 
suggest C(3) « 0.44 ,32^, slightly higher than the value predicted by the linear relation 
H54|l . The lower critical dimension for interfacial wandering induced by random-bond 
disorder is d^ = 5/3 |33l, below which they are unstable with respect to thermal 
fluctuations. 

The wandering exponent explicitly determines the fluctuation-dominated values 
of the critical exponents at the complete wetting transition. Complete wetting is 
pertinent to planar wall-vapour interfaces (say) with vanishing contact angle and 
refers to the divergence of the mean wetting layer (of liquid) thickness I, roughness ^j. 
and transverse correlation length ^|| as the bulk two-phase coexistence is approached. 
Writing the bulk ordering field h cx (fisatiT) — fj,) we define critical exponents according 
to the power laws 

I ~ , U - h-"^^" , C|| - /i-^r (55) 

The general theory of complete wetting is rather fully developed and can be understood 
using effective interfacial Hamiltonian models The usual starting point for the 

discussion of fluctuation regimes is the continuum model (O with a binding potential 

W{l)^hl + arP (56) 

where the Hamaker constant a is positive. Transfer matrix and renormalization group 
studies show that the critical behaviour falls in two scaling regimes '3^. The values of 
the critical exponents in these regimes can also be understood using a simple effective 
potential picture. As noted by Lipowsky and Fisher the form of the bending 
energy in the interfacial model suggests that fluctuation effects can be accounted for 
by minimising an effective potential 

WeffH) ^ hl + al-P + bl-^ (57) 




Figure 6. The vapour-liquid interface at complete wetting in a constrained 
geometry across the 2;— axis. See text for explanation. 



where t = 2{1/C{(j) — 1). The complete wettmg transition is therefore fluctuation- 
dominated, belonging to the so-called weak-fluctuation (WFL) regime, for r < p and 
leads to the critical exponent identiflcation 

with /3g° — v'^^ — C,{d)v'^° . This approach together with the flnal exponent 
identiflcation will be important in our discussion of filling in ordered and disordered 
systems. 



5. 2. From complete wetting to filling 

Complete wetting and 3D filling share a number of common properties. Both 
transitions exhibit two fiuctuation regimes reflecting the fact that there is always 
a relevant thermodynamic contribution to the appropriate effective potential. The 
difference in the values of the fluctuation repulsion exponents Tw and r can be traced 
directly to the extra lengthscale dependence in the effective stiffness at flUing. A 
relation between these exponents can therefore be found if we understand how an 
additional lengthscale in the stiffness alters the fluctuation-induced repulsion. To this 
end consider a complete wetting transition occurring in a flnite-size three-dimensional 
geometry with short-ranged forces (see Figure . The system is of infinite extent in 
the y direction (say) but of finite width L much bigger than the bulk correlation length. 
We may suppose that periodic boundary conditions apply across the system. In the 
presence of a flnite bulk field h the mean interface height is finite but diverges ash ^ Q 
with the growth of the film thickness falling into two regimes. If the width L is larger 
than the transverse correlation length ^|| ^ ft.^'^iif^) of a three-dimensional complete 
wetting layer of infinite extent, then the film grows with 3D-like critical exponents. 
If on the other hand L <C ^|| the growth of the wetting layer thickness must be 
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characterised by critical exponents pertinent to 2D complete wetting. Standard finite- 
size scaling ideas imply that the crossover between these regimes should emerge from 
the scaling hypothesis 

1{L, h) = /t-'3r(3)A(L;i''ii(3)) (59) 

where A(a;) is an appropriate scaling function. In the limit a; ^ oo the scaling function 
must approach a constant consistent with the bulk 3D limit. On the other hand as 
the argument vanishes we must impose 



A(a;)~a; "ii''' (60) 

in order to recover the correct 2D limit. Thus the asymptotic divergence of Z as ft- 
contains the finite-size dependence 

/3g°(3)-3g°(2) 

z^/i-^r(2)i (3) (61) 

This asymptotic divergence should also be understandable using a suitably modified 
effective potential measuring the free energy per unit length of the finite-size system. 
As /i ^ the interface behaves as a 2D-like complete wetting layer in an effective 
bulk field /le// oc hL and with modified stiffness '^eff oc SL. The divergence of / 
should therefore follow from the minimization of an effective 2D-like complete wetting 
potential 

Weff{h, L) = {hL)l + AL-"^/-^(2) (62) 

where the exponent (j) allows for the influence of the extra lengthscale dependence in 
the effective stiffness. In order to recover the correct finite-size result ()fil|l we must 
identify 

_ /3r(2) - /?r(3) _j_ _ 

/?-(2).||(3) +/3r(2) ^''^ 

This exponent is key to understanding the connection between interfacial 
wandering and 3D filling since it accounts for the correction to the standard fluctuation 
induced repulsion at wetting when the effective stiffness in the y— direction contains 
an additional length Se// = SL. This is also the case at filling where the effective 
stiffness along the wedge is proportional to the interfacial height 1^. Simple power 
counting implies that the fluctuation-induced repulsion exponent for the wedge is 

r„ = + t(2) (64) 

which reduces to 

In this way we can explicitly relate all the 3D wedge filling exponents in the FFL regime 
to the values of the interfacial wandering exponent for 2D and 3D planar interfaces. 
Thus Iw ^ ^ Cy™ with wedge wandering exponent 

^ (i + cS^-C(3)) ^^^^ 

which is valid provided both the two- and three-dimensional interfaces are well-defined 
and rough. That is the dimensionality must be such that d < du and d — 1 > d^- 
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Similarly, and subject to the same provisos, the critical exponent for the interlacial 
height reduces to 

- 2(1^ 

which is the central result of this section. The interpretation and implications of this 
result are discussed in detail below. 

5.3. Interpretation 

For pure systems the wandering exponent has values C(2) = 1/2 and C(3) = in two 
and three dimensions respectively. Consequently and recover the effective 
Hamiltonian predictions = 1/3 and — 1/4 discussed earlier. Further support for 
the connection between filling critical exponents and interfacial wandering arises when 
we generalise the above discussion to the d-dimensional wedge with d — 2 dimensions 
along the direction of translational invariance. The analysis is unchanged except that 
C(2) and ^(3) are replaced by ({d — 1) and respectively. Thus for a generalised 
wedge we predict 

- 2(1 -C(d)) ^ ' 
provided that fluctuations arc dominated by the d — 2 dimensional breather modes 
along the wedge. Substituting the thermal result ({d) = (3 — d)/2 recovers (|17|) for 
d < 3. For d > 3 (|t)8|l is no longer valid since the c?-dimensional interface is no longer 
rough. Notice that, in the limit c? — > 2, (|68|l recovers the correct two-dimensional 
filling result = 1 even though the breather-mode excitations are effectively zero- 
dimensional, and that C(l) = 1. 

For impure systems we are able to make some new predictions for filling in 
random-bond and random-field systems. For random-bond systems we use the values 
C(2) = 2/3 and C(3) = 0.44 to predict that in the FFL regime 

/3„ « 0.60 , w 0.60 , i^y w 1.1 (69) 

which are valid for sufficiently short-ranged forces. Assuming that there are only 
two fluctuation regimes and that there is a smooth cross-over from mean-field to 
fluctuation-dominated exponents we anticipate that the above predictions are valid 
provided the power law for the intermolecular forces p > 1.7. Thus dispersion-like 
(van der Waals) forces are irrelevant in the renormalization group sense and belong to 
the universality class of systems with short-ranged forces. Note that the value of the 
fluctuation repulsion exponent Tw is positive for both pure and random-bond systems. 

The case of 3D filling in systems with random-field disorder is less clear-cut. Using 
the values ({2) — 1 and ({3) — 2/3 we are led to the predictions Pw ~ = 3/2 and 
Uy = 2. However some caution is required in the assessing these predictions since the 
reliance on the marginal value C(2) = 1 is somewhat unsatisfactory. This situation 
is reminiscent of the two-dimensional limit of H68|) for filling with thermal disorder 
which the relies on the marginal value C(l) = 1. The fact that the breather-mode 
picture generates the correct two-dimensional result for this case leads some support 
to the above predictions for 3D random-field filling. However the situation is worse 
for random-field disorder since the fiuctuation repulsion exponent t^, is negative for 
d < 10/3. It is likely therefore that the exponent identification (|67|l is not appropriate 
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for 3D wedge filling with random fields. Here after we limit our discussion to thermal 
and random-bond disorder only. 

The above results for the values of the critical exponents for fluctuation- 
dominated three-dimensional wedge filling in pure and impure (random-bond) systems 
complement the known result for two-dimensional filling H20|l . It is therefore interesting 
to ask how the three-dimensional result (|67|l becomes the two-dimensional result as 
we compactify the wedge and reduce the number of dimensions along it. For pure 
systems with only thermal disorder the situation seems straightforward since the result 
H68|l is valid for 2 < d < 3. As noted above the generalised breather-mode result 
smoothly recovers the correct numerical result [3^ = 1 in this limit. There is however 
something remarkable about this since equating the two expressions (|^ and 
implies ^(2) = C(l)/2 yielding information about the allowed value of the wandering 
exponent. Noting that the lower critical dimension for pure systems is — 1 and 
that at this dimension C, = 1 we may conclude that equality of the breather-mode 
and wedge-covariant results for the critical exponents necessitates that for thermal 
fluctuations 



which is correct. Turning attention to random-bond systems however note that we are 
not permitted to extend the breather mode result down to d — 2 since the expression 
H68|l is not valid when d— l<di = 5/3 equivalent to d < 8/3. Therefore we are not 
able to equate distinct expressions for the filling critical exponents in two dimensions 
and are unable to discuss the values of these exponents for 2 < d < 8/3. However 
we are still free to discuss the properties of the filling critical exponents as d ^ 8/3 
from above. In particular, in view of the thermal result H7Q(I . it is natural to inquire 
what the value of the wandering exponent is for random-bond systems one dimension 
above the lower critical dimension. That is, what is the value of C(8/3)? As discussed 
earlier a reliable estimate for the value of the wandering exponent is provided by the 
linear relation H54|l which is only slightly below the exact result in two dimensions. At 
d = 8/3 this yields C(8/3) ~ 0.49 which we anticipate to be only slightly below the 
true result. This remarkable numerical coincidence leads us to conjecture that (|7(J|I is 
also valid for random-bond systems so that 



It is natural to ask whether this conjecture has any support from other examples of 
fluctuating surfaces. Lipowsky has considered interfacial wandering and unbinding 
within a generalised class of effective Hamiltonians with bending energy term of the 
forms JS] 



written in terms of the Fourier components of the collective co-ordinate. A remarkable 
property of such models is that for 77 < the fixed point behaviour (pertinent to 
short-ranged critical unbinding |36p is controlled by a single parameter rather than 
by r] and d separately. For our purposes we note only that the wandering exponent 
({d) = (3 — d — ?7)/2 and that the value at which this is unity allows us to identify 
dL = 1 — rj. This immediately implies that Cid-L + 1) = 1/2 consistent with the 
known result for interfaces in pure systems and our conjecture for interfaces subject 



CidL + 1) = ^ 



(70) 




(71) 




(72) 
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Figure 7. Schematic picture of the tilt fluctuations of the fiUed region in a section 
of a 3D wedge. 



to random-bond disorder. On the other hand the result l|7U() is not universaUy vahd 
since for interfaces in random fields ( ~ (5 — d)/3 implying C(c^l + 1) = 2/3. Further 
work is required to substantiate the conjecture (|71|) for random-bond disorder. 

6. Tilt and torsional mode fluctuations at wedge filling 

In this final section we return to the problem of wedge filling in pure systems and 
reconsider the transfer matrix analysis of the wedge Hamiltonian (jS)). Bednorz and 
Napiorkowski [37] have pointed out that some ambiguities in the construction of the 
infinitesimal transfer matrix arise because the bending energy term is proportional to 
the interfacial height. To avoid this they propose that the fluctuating field or order 
parameter appearing in the partition function functional integration is not lo{y) but 
rather z{y) = loiy)^^"^- The advantage of choosing this field is clear since the bending 
energy term is a simple Gaussian and one can map the problem immediately onto 
Euclidean quantum mechanics. Here we point out that this choice of order parameter 
emerges naturally when we consider the role of tilt and torsional fluctuations coupled 
to breather- mode excitations. 

Using the wedge Hamiltonian we wish to evaluate the partition function 



for a wedge of length L with flxed interfacial heights li = l{0) and I2 — 1{L) at the end 
points. The measure DIq denotes some suitable integration over the allowed values of 
the collective co-ordinate l{y) and will be discussed further below. The wedge model 
(O was constructed on the assumption that the only fluctuations that are relevant 
for filling correspond to local translations in the height of the filled region. Let us 
generalise this and allow for small, local tilts in the flat, filled region of the interface, 
that mimic the effect of the intrinsic capillary fluctuations around the constrained 
profile. That is we consider a restricted class of interfacial configurations, that within 
the filled region has the form 




(73) 



l{x,y) = lo{y) + (j)iy)x 



l{x,y) < a\x\ 



(74) 
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provided l{x,y) < a\x\. This is shown schematically in Figure The tilt angle is 
assumed to lie in some interval —ea<(j)<ea with < e < 1. Outside of the filled 
region we assume that the interface remains within a microscopic distance from the 
wall so that l{x,y) » a\x\. Substitution into the capillary- wave model leads to a wedge 
Hamiltonian that allows for breather, tilt and torsional modes and is a functional of 
the two fields Iq and (p. It is straightforward to show that the breather-tilt-torsional 
(BTT) model has the form 

where the generalised wedge potential is 

VilQ,^)^^{0^ + (f>^)lo + Cll~P (76) 

The bending energy coefficient for the breather mode Ki = ^ and is unchanged 
from the simpler model ©• The coefficients appearing in the torsional terms are 
K2 (x Sa~^ and (x 'Sa~^ respectively. The precise values of these coefficients 
are not important in the critical regime. It is also convenient to re-write the wedge 
potential as 

V{lo,cl>)=tlo + Cll-P + -^Ho (77) 

a 

Notice that when we set the field (/)(?/) = at all points along the wedge the BTT 
Hamiltonian reduces to the wedge model (|85|l . To assess the relevance of tilt and 
torsional modes on the critical behaviour of the model we consider a rescaling of the 
co-ordinates and fields similar to that discussed earlier. With the additional tilt field 
we now consider the properties of the model under the mapping 

y^y' = l , lo^l'o = ^ , <j,^4>' ^^h-'i- (78) 

where, as before, the wedge wandering exponent C, — 1/3. Under this rescaling the 
coefficient of the breather- mode bending coefficient Ki and the term (jy^l^ appearing 
in the wedge potential remain invariant. The linear temperature-like scaling field t 
and Hamaker constant C rescale as earlier (see (|12|l ') whilst the new tilt and torsional 
bending energy coefficients renormalise to 

6-4/3^2 , K'^^b-^'^K^i (79) 

and are therefore irrelevant for all ranges of intermolecular force. This means that 
they do not affect the critical behaviour and that the BTT model can be simplified to 

HBTT[k,<f\= + ^^^^ 

The tilt field (j) is therefore non-interacting but still coupled to the breather-mode 
excitations. As the interface unbinds from the wedge bottom the tilt fiuctuations 
become increasingly massive and are suppressed. An important consequence of this 
is that when we integrate out the tilt field the particular choice of the parameter e 
becomes unimportant in the scaling limit. 
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The final ingredient in the construction of the model is the nature of the 
measure Dl^Dcj) appearing in the functional integration over the Boltzmann weight 
^-f)HBTT[iM ^ where (3 = l/ksT. In reducing the dimensionality of the collective co- 
ordinate l(x,y) {l{y) , 4>{y)) some care is required in defining the new measure 
since we must impose the correct normalization conditions. This is most easily 
illustrated for 2D wedge filling in which there is no y-direction and the choice of 
measure DIqDcI) relates only to one / and one 4> variable. For this case the Boltzmann 
weight for a tilt configuration is simply e"''^^'''^-' and determines the un-normalised 
probability P(/oj 4') for finding an interface with tilt angle (j) and mid-point height Zq. 
Integrating over the allowed values of determines the mid-point height PDF P{Iq). 
The properties of this PDF are known from exact calculations. In particular for short- 
ranged forces (belonging to the filling fiuctuation regime) the function has the simple 
exponential form P{Iq) oc g-'^P^i^-'^)'-^ ^j^g scaling limit. Thus we require that as 

/ea 
i^/od0e-^^(^'-"'+^')'°/" (81) 

which necessitates that we interpret the 2D measure DIq cx ^/TodlQ. In other words 
the correctly normalised PDF P{lo,(jji) oc y/loe^^^^-'^K An equivalent argument is 
that after integrating over the tilt fluctuations the resulting PDF P{lo) must have the 
correct short-distance expansion as the interface approaches the wall. For 2D systems 
with short-ranged forces the PDF has the short-distance expansion P{Iq) ~ with 
7 — ^ 1- Recall that for pure systems (3^ — ^ (in 2D) implying that the two- 

field PDF P{Iq,4>) must contain a ^/Tq prefactor in-order to preserve thermodynamic 
consistency. This argument is supported by explicit results for the 2D wedge filling 
for contact binding potentials (see [Appendix A| ). 

A similar line of reasoning applies in higher dimensions. To determine the 
form of the 3D measure we first focus on the filling transition occurring in a cone 
geometry. The reason for this is that the breather-mode excitations are essentially 
zero-dimensional and we only have one height variable Ic measuring the height of the 
interface above the cone vertex. Up to a normalization factor the probability of finding 
the interface at height Iq with tilt (j) is determined by the Boltzmann weight e~^^'' 
where Vc is the free-energy cost of the configuration. The leading order term in the 
cone binding potential contains a term proportional to the area 1^ of the filled region. 
The prefactor of this is f3T,{d'^ — -I- 0^) similar to the 2D wedge and integrating over 
the tilt modes (and a trivial rotational degree of freedom) determines the mid-point 
height PDF Pc{lc)- In order to preserve the simple Gaussian form of this function 
we must interpret the 3D measure Die oc Icdlc- For the 3D wedge we discretise the 
system in the ?/— direction and replace the fields {l{y)} and {<j){y)} by sets of continuous 
variables Ij and 0j. The wedge measure must be proportional to a product of cone 
measures and thus we interpret the functional integral for the BTT Hamiltonian as 

jDloD(t>^ J J J Y[lo{j)dloij)dcPj (82) 

On integrating over each tilt angle we recover the wedge model © and identify the 
functional integral appearing in (|73|l as the continuum limit of 

/ /••■/n^o(j)''''^'oO') (83) 
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In terms of the variable Zj — loiiY^"^ the functional integral therefore becomes 

\{dz{j) (84) 

j 

and has a well-defined infinitesmal transfer matrix limit. In terms of the new variable 
the continuum wedge Hamiltonian is 

where V{z) = V{z'^/'^ The rest of the analysis is identical to that described by 
Bednorz and Napiorkowski and the problem can be mapped immediately on to 
quantum mechanics. Thus the partition function may be written as the spectral 
expansion 

ZihM-.L) = MhYMWe-"''"'' (86) 

n 

where the eigenfunctions and eigenvalues satisfy the Schrodinger equation 
9a 



IQP2^M^)" + V{z)Mz) = EMz) (87) 

with boundary conditions "0(0) = ^j{oo) = 0. Hence the PDF for finding the interface 
at height Iq can be identified with Pw^Jo) — |V'o(^^^^)P whilst the wedge free energy 
and correlation length follow from — Eq and — (3{Ei — Eq) respectively. For 
the choice of wedge potential IjlUI) it is straightforward to deduce the existence of two 
fluctuation regimes with the critical exponents discussed earlier. 

Our final remarks concern the form of the interfacial height PDF P(/o) in the 
FFL regime. The intermolecular forces are irrelevant for p > A and in the limit 9 a 
the scaling properties of the PDF follow from solution of the differential equation 

- Y^MzY' + 2S(0 - a)ziMz) = EoM^) (88) 

It follows immediately that the wave-function and hence PDF are scaling functions of 
the dimensionless variable u — y/]JS{9/a — IY^^Iq and we can write 

/n \i/4 

P(?o) = V^f^-lj g{^{9/a-iy/^lo) (89) 

where g{u) is a suitably normalised scaling function such that J dug{u) = 1. The 
scaling of the PDF is, of course, to be anticipated since in the FFL both l^ and 
diverge with the same critical exponent. We are not aware of a closed form solution 
to the differential equation in terms of elementary functions. However both the 
asymptotic short-distance and large-distance behaviours agree with thermodynamic 
and sum-rule requirements. First consider the large distance limit corresponding to 
u ^ 1. From solution of (|88|) it is straightforward to show that the PDF decays like 



P{h) ^ ^{9 /a - l)l/4e-2;3SV2(e/a-l)i^ (qq) 



up to terms of order in the argument. This is precisely in keeping with macroscopic 
requirements. When the interface is far from the wall the PDF can be identified with 
g-/3AF ■y^f]^g]-g /\^p jg ^j^g free-energy cost of a macroscopic configuration constrained to 
pass through the height Zq- In the wedge geometry this macroscopic configuration has 
the shape shown in FigureQand the free-energy cost is dominated by a thermodynamic 
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Figure 8. Rescaled probability distribution function for the interfacial height 
above the wedge bottom found from the breather mode model. Integration of the 
PDF determines the density profile in the wedge. 



area-like contribution. This can be readily calculated making use of Young's equation 
yielding AF = 2T,^/2{e - a)/all. The short-distance expansion (SDE) of the scaling 
function g{u) corresponding to the limit w ^ 1 is also consistent with known critical 
exponent relations. For 2D and 3D wedge filling the scaling function for the PDF has 
the short-distance behaviour 

g^u-* (91) 

where the SDE critical exponent 7 = I//?™ — 1. This exponent relation follows from 
rather general sum-rule and scaling arguments. Thus for 3D filling in pure systems 
with /3m = 1/4 we anticipate 

P{h) - - 1) I'o (92) 

for \/ j3Ti{9/a — IY^^Iq <^ 1. This emerges naturally from solution to the Schrodinger 
equation since in terms of the z variable ipi^) ~ -z as z — > and recall P{lo) — 
IV'('o^^)P- The large-distance and SDE properties of the PDF give strong support to 
the correctness of the one-dimensional wedge Hamiltonian model. The numerically 
determined scaling form of the PDF found from solution to the Schrodinger equation 
is shown in Figure |S1 

7. Conclusions 

In this paper we have addressed a number of issues pertaining to 3D wedge filling 
transitions. We began with a discussion of wedge covariance relations between filling 
and critical wetting and showed that the fluctuation-induced covariance observed for 
2D filling has a mean-field analogue. Similar to the 2D case classical covariance 
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appears to be specific to systems for which the critical wetting specific heat exponent 
as = 0. Numerical studies of mean-field filling within a Landau theory show that 
classical covariance is not limited to the shallow wedge limit. Importantly this 
provides a constraint on the possible form of effective Hamiltonian models of wedge 
filling and implies that a model suggested by Rejmer et al [7] is not consistent with 
more microscopic approaches. It is likely that in order capture the correct classical 
covariance for wetting and filling it will be necessary to construct non-local interfacial 
models. This will be discussed elsewhere. 

In the second part of our paper we turned attention to fluctuation effects 
at 3D wedge filling. Our central new result is that the critical exponents which 
characterise the fluctuation-dominated regime can be related to the values of the 
standard wandering exponent for planar interfaces in 2D and 3D (bulk) systems. 
The expressions for f3w and the wedge wandering exponent C^^ expressed in terms of 
C(2) and C(3) are rather general and allow the discussion of filling in systems with 
random-bond disorder. A surprising feature emerging from our study concerns some 
intriguing numerical coincidences for the value of the bulk wandering exponent C,{d) 
as one approaches the lower marginal dimension for wedge filling in pure and random- 
bond systems. In particular we are led to conjecture that for random-bond disorder 
C(8/3) = l/2. 

In the final part of our paper we addressed some queries that have been highlighted 
concerning the pseudo-one-dimensional wedge Hamiltonian theory of 3D filling in pure 
systems. Bednorz and Napiorkowski have suggested that these can be avoided 
if one supposes that the fluctuating field is not the interfacial height lo{y) above the 
wedge bottom but loijjY^^. We showed that this choice indeed emerges naturally 
if one considers the coupling between tilt and torsional fluctuations described by 
a second field and breather-mode excitations of l^iy). The form of the PDF 
for the interfacial height predicted by the wedge Hamiltonian model is shown to 
have the correct short-distance and large-distance behaviour dictated by macroscopic 
arguments and critical exponent relations. We believe this strongly supports the 
internal consistency of the effective wedge Hamiltonian model. As mentioned in the 
introduction the values of the critical exponents predicted by the wedge Hamiltonian 
are in very good agreement with Ising model simulation studies. A more stringent test 
of the theory would be to study the scaling of the interfacial height PDF. This can 
be readily extracted from Ising model studies since it corresponds to the derivative of 
the magnetization profile measured along the vertical above the wedge bottom. To do 
this quantitatively it would be necessary to generalise the present theory to study the 
finite-size scaling of the PDF and magnetization profile in the same geometry studied 
by Milchev et al ^1 . This will be a topic of future work. 

Finally we mention that ideally one would hope to circumvent the entire wedge 
Hamiltonian theory and replace it with a more microscopic approach. In particular it 
would be extremely useful if one could apply a renormalization group analysis to derive 
the wedge filling critical exponents from the full capillary- wave model (U)). To do this 
one would have to identify a new fixed-point distinct from the Gaussian (or strong- 
fiuctuation regime) fixed-point used to study wetting transitions. Such additional fixed 
points for the capillary-wave model may exist if one generalises the spatial rescaling 
to account for the anisotropy of correlation lengths across and along the wedge. One 
would hope, of course, that such an approach would lead to the same predictions 
for universal critical exponents and scaling functions found using the effective one- 
dimensional wedge Hamiltonian. As well as being an independent check on the validity 
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of the wedge Hamiltonian model (in the scahng hmit) such an approach would also 
shed much light on the emergence of an effective one-dimensional theory from a higher 
dimensional model. 
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Appendix A. Tilt mode fluctuations at 2D wedge filling 

In this Appendix we provide further justification for the role played by tilt fluctuations 
at 2D wedge filling making contact with exact results obtained from interfacial 
Hamiltonian studies. We restrict attention to systems with short-ranged forces 
corresponding to the FFL regime and set kgT — 1 for convenience. Consider 
the 3— point probability distribution function P^{li,—x;lo,0;l2,x) determining the 
probability for finding the interface at heights h, Iq and I2 at positions —x, and x, 
respectively. Exact transfer-matrix results for a continuum interfacial model of filling 
in shallow wedges show that Pwih, —x; Iq, 0; I2, x) decomposes into the product 38 



where PwHo) — 21]{d — a)exp[— 2S(6' — a)lo] is the wedge midpoint 1— point PDF, 
I — I — ax is the local relative height to the wall (recall that x is positive) and 
P^{12,X2\IitXo) is the planar conditional probability distribution function I38| : 



If the wedge midpoint height is much greater than the mean wetting height for a 
planar substrate, i.e. Iq 3> 1/2I]0, which is fulfilled in the asymptotic critical regime 
and the position x < lo/9, the second term in Equation ljA.2p can be neglected and 
the conditional probability distribution function becomes a gaussian distribution of 
average Iq — Ox and dispersion <t^ = x/T,. The interface will be affected by the presence 
of the substrate when Iq — Ox^ y^x/T,. This expression defines the effective value of x* 
at which the interface is in contact with the substrate. We consider the angle formed 
by the horizontal and the line that joins the interface position at the midpoint and 
the position (relative to the substrate) at x* (see Figure IXT|) . The tilt angles A^i and 
A62 are defined as the departures of the previously defined angles with respect to the 
average value 9. We stress that the interface is not considered stiff but is allowed to 
wander due to capillary-wave fluctuations. However, the length scale on which these 
fluctuations occur is much smaller than the one that characterizes the constrained 



PUk,~x;lo,0;l2,x) = PUlo)Pi;{lux\lo,0)P^{l2,x\lo,0) 



(Al) 




(A.2) 



profile. When A6i,2 < 6* < 1, then li,2 - Iq + Ox* 



x*M 



1,2- 




Figure Al. Plot of the typical relative intcrfacial position 1 = 1 — a\x\ close to 
the 2D wedge filling transition. A^i and A02 define the tilt fluctuations of the 
interface at x* . See text for explanation. 



Under these assumptions, the PDF for a wedge midpoint position and tih 
angles A9i and A6'2 can be approximated as 

P(/o, « 2E(0-a)e-2S(''-")'°— e-'^[(^''^)'+(^''^)'](A.3) 

2'K 

where the tilt angle ranges can be extended safely to the complete real axis. We define 
the tilt angles (j) and (j)' as (A^i + A02)/2 and {AOi — A02)/2, respectively. Taking into 
account that in the scaling limit x* « Iq/^ ~ ^o/o^j Equation HA.3() can be expressed 
as: 

X J e = (A.4) 

V na 

Integrating all the values of 0', we obtain: 



«^^!^(0-a)e-^(«^-"^+^^) (A.5) 

y'TTa 

in agreement with Equation H81|) and with the correct measure DIq ^ ^/IqcIIq. 
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